The activation cascade of the broad-spectrum antiviral bemnifosbuvir characterized at atomic resolution

Bemnifosbuvir (AT-527) and AT-752 are guanosine analogues currently in clinical trials against several RNA viruses. Here, we show that these drugs require a minimal set of 5 cellular enzymes for activation to their common 5′-triphosphate AT-9010, with an obligate order of reactions. AT-9010 selectively inhibits essential viral enzymes, accounting for antiviral potency. Functional and structural data at atomic resolution decipher N6-purine deamination compatible with its metabolic activation. Crystal structures of human histidine triad nucleotide binding protein 1, adenosine deaminase-like protein 1, guanylate kinase 1, and nucleoside diphosphate kinase at 2.09, 2.44, 1.76, and 1.9 Å resolution, respectively, with cognate precursors of AT-9010 illuminate the activation pathway from the orally available bemnifosbuvir to AT-9010, pointing to key drug–protein contacts along the activation pathway. Our work provides a framework to integrate the design of antiviral nucleotide analogues, confronting requirements and constraints associated with activation enzymes along the 5′-triphosphate assembly line.


Introduction
The recent Coronavirus Disease 2019 (COVIDAU : Pleasenotethat}COVID À 19}hasbeenfullyspelledo -19) crisis due to Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has evidenced the need for safe and potent antivirals, in conjunction with accurate and rapid diagnostics.First, at the public level, prophylaxis around index cases may significantly curb emergence by cutting viral load and transmissibility.Second, at the individual level, early therapy may limit virus spread to vital organs, avoiding complications.
Nucleoside/nucleotide analogues (NAs) represent the first family of compounds that have been used as direct-acting antivirals several decades ago [1][2][3] against DNA viruses such as herpesvirus, hepatitis B virus (HBV), the retrovirus human immunodeficiency virus type 1 (HIV-1), and RNA viruses such as hepatitis C virus (HCV) and SARS-CoV-2.Upon reaching their target cells, antiviral drugs often act after intracellular activation.Early biologically active nucleoside analogues (NAs) such as ara-C (cytarabine) penetrate the target cell efficiently and are transformed into the nucleoside analogue 5 0 -triphosphate (NA-TP) by a series of kinase reactions (reviewed in [2][3][4]).
ProTide prodrug technology has been implemented to deliver NA carrying a 5 0 -phosphate or phosphonate [5].Aryloxy phosphoramidate prodrugs of, e.g., tenofovir alafenamide (TAF) for HIV-1 and HBV, sofosbuvir (SOF) for HCV, and remdesivir (RDV) for SARS-CoV-2 have met clinical success [5][6][7].They achieve cell penetration through shielding the phosphate charges, and preferential hydrolysis at the P-N bond bypasses the first-often limiting-nucleoside kinase reaction, allowing significant building up of NA-TP pools.NA-TPs poison viral RNA synthesis specifically, provoking premature viral RNA chain termination or chemical/genetic corruption of the viral nucleic acid [8].
Two main parameters determine NA potency: the concentration ratio of NA-TP over its natural NTP competitor, and the selectivity of the viral RNA-dependent RNA polymerase (RdRp) for use of the NA-TP [9].NA-TP pools must be formed efficiently upon intracellular activation, and the differences of the NA-TP scaffold relative to its natural NTP counterpart must remain "below the radar" of the viral RdRp.
General knowledge, structural and functional data about cellular enzymes along a given NA activation pathway are fragmented.Numerous NAs have been designed showing appropriate poisoning of viral RNA synthesis through their 5 0 -TP in vitro, but clinical development has failed for lack of transmembrane permeability, intracellular metabolic activation, and/or cellular toxicity.Clearly, NA drug-design needs an integrated view (i.e., structural, functional, and mechanistic) from the delivered NA up to the ultimate inhibited viral enzyme reaction accounting for antiviral effect, as well as the molecular basis for (non)interaction with cellular enzymes.
Such a global picture is emerging for 2 FDA-approved NA prodrugs directed against RNA viruses, SOF and RDV (S1 Fig) .The activation pathway of RDV has been elucidated and profiled in a variety of tissues [10].RDV is an aryloxy phosphoramidate prodrug (GS-5734, Veklury) converted in 2 steps to the 5 0 -monophosphate NA by esterases belonging to 2 families (Cathepsin A (CatA) and/or carboxylesterase 1 (CES1), followed by the cleavage of the P-N bond by a histidine triad nucleotide (HINT) phosphoramidase [11].The resulting monophosphate NA is then converted to RDV-TP by the subsequent action of 2 cellular phosphotransferases.However, although detailed studies at the atomic level exist to understand accommodation of the 1 0 -CN group at the SARS-CoV-2 RdRp active site [12,13], structural insight of the relevant activation intermediates is lacking for all enzymes in the activation pathway.
SOF (PSI-7851, Sovaldi) undergoes the same CatA/CES1 deprotection pathway as RDV but is activated to the 5 0 -triphosphate by the UMP-CMP kinase and the nucleoside diphosphate kinase (NDPK) sequentially [14].Here again, although the mechanism of chain-termination at the HCV NS5b RdRp has been elucidated at atomic resolution [15], the structural basis of activation remains uncharacterized.
Bemnifosbuvir (AT-527 (hemisulfate salt), AT-511 (free base)) is a guanosine analogue currently in clinical trials against SARS-CoV-2 and HCV.It showed in Phase II clinical trials a 71% risk reduction in outpatients with moderate COVID-19 (MORNINGSKY; NCT04396106) [16,17], but the associated Phase III study ended prematurely as it did not meet its primary end point [17].It is currently under investigation in a global Phase III clinical trial in outpatients at high risk for disease progression (SUNRISE-3; NCT05629962), as well as evaluated as an anti-HCV drug [18] in combination with the NS5A inhibitor ruzasvir (NCT05904470).Its epimer AT-752 is currently in clinical Phase II against Dengue virus (NCT05466240) [19,20].
Most information relative to activation pathways has been obtained through measurements of intracellular concentrations of NA intermediates [10,11,14,24].Studies described herein aim to ascertain which individual enzymes are involved in the activation pathway, as well as clarify their structural and functional mode of interaction with activation intermediates.Interaction maps of NAs with both activating partners and viral enzymes at atomic resolution would guide improved design of novel NAs.
In this work, we have identified a minimal set of 5 human enzymes involved in the activation of bemnifosbuvir and AT-752 to their active 5 0 -TP form AT-9010. We demonstrate distinct stereo-preference of CatA and CES1 for AT-527 and its epimer AT-752, respectively.We elucidate an ordered activation pathway with potential limiting steps, including substrate specificity requirements at the N 6 -purine position for ADALP1.Structural analysis with key human enzymes HINT1, ADALP1, GUK1, and NDPK in complex with the bemnifosbuvir/ AT-752 intermediate NA at 2.09, 2.44, 1.76, and 1.9 Å resolution informs on both the N 6 purine and 2 0 -ribose chemosteric requirements leading to active NAs targeting viral RdRps.

Results
The putative metabolic pathway of AT-511 (free base form of AT-527, bemnifosbuvir) and AT-281 (free base form of AT-752) is shown in Fig 1A.
CatA and/or CES1 presumably precede a spontaneous decomposition yielding the L-alanine phosphoramidate prodrug AT-551, common to both AT-752 and AT-527 pathways.AT-551 is then subjected to 4 enzyme-mediated reactions leading to the active 5 0 -triphosphate form, AT-9010.Apart from the N 6 -deamination reaction, this pathway has been inferred from that of the anti-HCV drug SOF, which carries a uracil base instead of the N 6 -methyl diamino purine base (a guanine precursor) of either AT-752 or AT-527.The SOF pathway, however, has been determined indirectly, by measuring intermediates in cells treated with SOF [11,14].

Five enzymes are needed to convert AT-527 (bemnifosbuvir) or AT-752 to their active form AT-9010
We expressed, purified, and crystallized 4 of the human enzymes widely distributed in human tissues (https://proteinatlas.org)and supposedly involved in the prodrug activation pathway (Figs 1A and S2).Together with human enzymes CatA and CES1, we challenged HINT1, ADALP1, GUK1, and NDPK with their anticipated substrates depicted in Figs 1A and S3.We made use of an HPLC-UV method to follow the conversion of the substrate, together with a kinetic monitoring relative to known standards (Fig 1B).
Both AT-527 and AT-752 can be enzymatically converted to AT-9010 by an ordered series of reaction involving nonspecific esterases CatA/CES1, HINT1, ADALP1, GUK1, and NDPK (Table 1).For AT-511, nonspecific esterases CatA and CES1 exhibit an approximatelyAU : Pleasenoteth 300-fold difference in activity, whereas hydrolysis efficiency is quite similar for AT-281 (<3-fold difference).Both for control and AT compounds, the early steps of the pathway are the slowest, with di-and triphosphate formation showing the fastest turnovers.

The phosphate stereo-selectivity of CatA/CES1 usage
The phosphorus atom of the aryloxy phosphoramidate moiety is chiral, resulting in 2 possible epimers.AT-511 (S P isomer) is the free base form of AT-527, and AT-281 (R P isomer) is the free base form of AT-752 (S3 Fig) .As shown for SOF, the difference of antiviral activity between the 2 stereoisomers could be due to the activation of the prodrug and more specifically to the stereo-selectivity of the enzymes involved in the first hydrolytic step.
CatA or CES1 were incubated with either AT-511 or AT-281, and comparative velocities of AT-551 formation were measured (S4 Fig) .CatA always shows higher activity than CES1.CatA prefers by 14-fold the S P isomer AT-511 over the R P isomer AT-281 (9.8 ± 0.3 versusAU : Pleaseno 0.7 ± 0.1 nmol/min/nmol protein).Interestingly, however, a reversed stereo-selectivity is observed for CES1 relative to CatA: CES1 hydrolyses AT-281 approximately 10-fold faster than AT-511 (0.27 ± 0.02 versus 0.033 ± 0.004 nmol/min/nmol protein, respectively).These results align with those obtained with a variety of clinically relevant prodrugs; for SOF, TAF, and RDV, CES preferentially hydrolyzes the R P isomer, whereas CatA prefers the S P isomer of the aryloxy phosphoramidate moiety [10,11,14].We conclude that nonspecific esterases of the CatA/CES1 type are able to convert bona fide McGuigan ProTides into HINT1 substrates from both AT-752 and AT-527.The respective abundance of each enzyme in cells or tissues might be useful as a predictor of ester hydrolysis efficiency in these given cells or tissues.

Specificity and order of the reactions
AT-551 could not act as a substrate for ADALP1 (S5 Fig) .This result indicates that the L-alanine phosphoramidate precludes efficient either binding or catalytic conversion of AT-551 with ADALP1 (see structural data below).AT-8003 could neither be a substrate for GUK1, indicating that the N 6 position must be either unsubstituted or only an O 6 to allow binding and phosphorylation of the NMP (S5 Fig) .Unless other enzymes exist that could use these substrates and bypass the proposed activation sequence, the absence of reactivity reported here is consistent with the order of the reaction pathway reported in Fig 1A .Structural and functional analysis of HINT1 substrate specificity HINT1 belongs to the HIT (histidine triad) protein superfamily and has been shown to hydrolyse the P-N bonds in SOF and RDV [10,11,14].HINT1 is able to hydrolyze AMP-NH 2 (control) into AMP, and AT-551 into AT-8003 with specific activities of 101 ± 8 and 0.7 ± 0.2 nmol/min/nmol protein, respectively (Fig 1B and Table 1).HINT1 is thus approximately 140-fold faster with the control compound than with AT-551.
For AT-551, reaction occurred readily under either crystal soaking or co-crystallization conditions.Co-crystallization resulted in AT-8003 bound to the active site on 1 monomer (Fig 2A).The crystal structure (2.09 Å resolution) shows a homodimer, each monomer being almost perfectly superimposable (RMSD approximately 0.18 Å on the visible 126 aa).The bound product pose is significantly different from all existing HINT1 structures [21] in complex with AMP, GMP, or NMP-phosphoramidate ligands (Fig 2B).In all available structures so far, the NA is buried in a pocket, contacting conserved amino acids Ile 44 and Asp 43 involved in purine base stacking and ribose 2 0 ,3 0 -OH ribose bi-dentate binding, respectively.Here, the methylated diamino purine base is bound outside the nucleobase binding pocket, shifted by 90˚and making a single hydrogen bond with Ser 45 that could occur with any purine or pyrimidine base.Trp 123 is facing the ribose offering a weak hydrophobic stabilization.Remarkably, the 2 parts of the AT-8003, which differentiate it from guanosine, are not engaged in any amino acid contact.Both the N 6 -amino purine and the 2 0 -C-methyl, 2 0 -F are projecting into the solvent.
The phosphate is closely superimposable to that observed in other published structures [21].The catalytic His 112 nucleophilic nitrogen is 3.1 Å away from the phosphate, and the hydrolytic water molecule is positioned opposite to the putative leaving amine and general acid His 114 and/or His 51.In a similar complex (PDB: 5IPE), the O 6 -guanine of GMP is located at 3.3 Å from the Ile 18 side chain, and there is space to accommodate an N 6 -methyl amino group at the solvent-protein interface.The AT-8003 ribose pucker is of Northern (N) configuration: C3 0 up and C2 0 down (the Southern (S) pucker has the C3 0 down and C2 0 up), as this conformation must be imposed by the presence of its 2 0 -C-methyl.The ribose pucker of GMP and AMP analogues in all other structures being S, we surmise the 2 0 -C-methyl is driving this alternate preferential NA binding mode seen in all structures reported here.
Under our experimental conditions, ADALP1 converts only nucleoside monophosphate substrates: No activity was detected with the nucleosides AT-229 and AT-259.positions of AMP (light green-PDB: 5ED6), CMP (orange-PDB: 5KM2), UMP (yellow-PDB: 5KM3), GMP (purple-PDB: 5KM1) in respect to AT-8003 (white) in the HINT1 active site.AT-8003 is positioned almost orthogonally to other nucleotides bound in other structural models, and only the phosphate positions are undistinguishable.The right panel is the same figure turned 45˚.(D) Human ADALP1 analysis and specificity.Activity of ADALP1 with N 6 -Me-AMP as control substrate, AT-8003 as anticipated substrate and N 6 -modified analogues of AT-8003: AT-8004 with a N 6 -dimethyl group, AT-8010 with a N 6 -npropyl group and AT-8002 with a N 6 -amino group.About 100 nM ADALP1 (or 1 μM for AT-8002) was incubated 2 h at 37˚C with 200 μM substrate.Bars show mean values (± SD) of 3 independent experiments.(E) ADALP1 in complex with AT-8001.Human ADALP1 structure is represented in ribbon and transparent surface, while compound is represented in sticks with heteroatom colors.The molecule is deeply buried in a closed conformation of ADALP1.On the right side, a zoomed part of the molecule as well as the zinc ion (gray ball) is presented in the active site.Side chains of interacting residues with the molecule are shown in sticks.The F o -F c omit map of the molecule is presented at 1 σ.The data underlying this figure can be found at https:// doi.org/10.2210/pdb8pwk/pdband https://doi.org/10.2210/pdb8qch/pdb.https://doi.org/10.1371/journal.pbio.3002743.g002 Interestingly, the reaction is approximately 45-fold slower with AT-8002 as substrate compared to the natural substrate N 6 -Me-AMP, indicating that at least one substituent on the nitrogen atom may enhance the leaving group ability.
We co-crystallized ADALP1 and AT-8003.The solved ADALP1 structure adopts an α/βbarrel architecture similar to that of ADAs [29], with the essential Zn 2+ ion tightly bound in the vicinity of the catalytic center (Figs 2E and S6A).Fourteen α-helices surrounding 8 parallel β-strands constitute a TIM-barrel fold [30].The asymmetric unit contains 8 copies of the protein.Since ADALP1 was purified as a homogeneous monomer (S6B Fig), the octamer likely appeared during crystallogenesis.Although ADALP1 was co-crystallized with the substrate AT-8003, the product AT-8001 was found bound at the active site.The compound is found in all 8 chains, and all the chains found in the asymmetric unit closely superimpose (RMSD of approximately 0.26 Å) except for the minor loop and N-and C-terminal ends (S6C Fig) .The structural zinc is coordinated by side chains of 4 residues (His 24, His 26, His 208, and Asp 293) and a water molecule, in a trigonal bipyramidal geometry [31].Water bridges His 232 and the resulting O 6 of the purine base, suggesting that both His 293 and the water molecule did play a role in the hydrolysis of the amine.AT-8001 is buried in a closed hydrophobic pocket that can only accommodate a monophosphate compound checked and stabilized by Asn 28, His 74, Ser 106, and Thr 107.
Unlike in the HINT1:AT-8003 complex, there is no ambiguity in substrate nor product binding pose [30].Once the substrate is bound, the catalytic water attacks the purine C 6 , which likely stays at this location while the methylamine is leaving in the opposite exit channel.

Structural and functional analysis of GUK1 specificity
GUK1 is a nucleoside monophosphate kinase able to transfer a phosphate group from ATP to a guanosine 5 0 -monophosphate, yielding GDP and ADP [25].Simultaneous binding of monophosphate substrate and triphosphate donor precedes conformational closure and phosphate transfer.This enzyme metabolizes several purine analogues such as abacavir and ganciclovir [32,33].Expectedly, GUK1 is able to phosphorylate GMP into GDP, but also AT-8001 into AT-8500 (Fig 1B and Table 1).GUK1 is much faster with GMP (control) than with AT-8001 (respectively, 6,800 ± 300 nmol/min/nmol protein versus 34 ± 5 nmol/min/nmol protein).GUK1 exhibits the highest turnover in the activation chain (Table 1), but the 2 0 -C-methyl, 2 0 -F modification slows it down approximately 200-fold.
GUK1 was crystallized and its structure solved at 1.76 Å resolution in complex with AT-8001 (Fig 3A and Table 2).
The overall fold of GUK 1 closely resembles that of the yeast and mouse guanylate kinases [34,35].Two β-sheets and 8 α-helices, the last 2 separated by a π-helix, define 3 structurally and functionally distinct subdomains: the core, the NMP-binding site, and the lid [35].The core is defined by residues 5 toAU : Pleasenotethaten À dasheshavebeenchangedtotoforclarityinnumber 31, 97 to 123, and 165 to 194 (helices α1, α4, α7, and α8; strands β1, β7, β8, and β9).The NMP-binding site is a mobile structural element defined by residues 37 to 89 (starting with a protruding loop and followed by helices α2 and α3; strands β3, β4, β5, and β6).The lid where the phosphate donor (e.g., ATP) binds is also mobile and defined by residues 126 to 156 (helices α5 and α6).These parts are interconnected with 4 hinges that allow conformational changes from an open to closed state [34].The structure presented here is an open conformation (Fig 3A ), with the AT-8001 substrate binding through extensive contact with both the phosphate moiety and the nucleobase.Thus, on the one hand, Tyr 81, Tyr 53, and Arg 41 share hydrogen bonding with the oxygen of the phosphate, while Arg 44 makes a contact through a water molecule.On the other hand, Glu 72 and Ser 37 engage hydrogen bonds with the base, while Tyr 53 contributes to positioning the base through π-stacking.The protein structurally checks the substrate by contacting both the purine base and the presence of a monophosphate, with no significant check on the ribose at this stage.We modeled AT-8001 into a closed GUK1 conformation based on an energy-minimized model of hGUK1 (S7 Fig) .The base is constrained by interactions with Ser 37, Glu 72, and Thr 83.The major difference observed with bound NMP substrates is the ribose pucker conformation.Again, in the structure with bound GMP, the 3'-OH ribose is in S pucker position to accommodate the closed conformation of the active site, but not for AT-8001, as described for HINT1 above.To avoid a steric conflict of the 2 0 -C-methyl, the compound has to slide into the core where the base is π-stacked by Tyr 81 and stabilized by Ser 37 and Glu 72.Both Arg 41 and Arg 44 contact the phosphate; Arg 137 and Arg 148 are expected to catalyze phosphate transfer from ATP bound in the lid, but the repositioning creates distances not compatible for optimal transfer.This positional adjustment likely explains the approximately 200-fold discrimination of AT-8001 versus the natural substrate GMP (Table 1).

Structural and functional basis of NDPK specificity
NDPK is able to transfer a phosphate group from a triphosphate nucleoside donor to a nucleoside diphosphate acceptor [26,36], to yield the active triphosphate form of bemnifosbuvir AT-9010.The reaction proceeds through a so-called "ping-pong" mechanism in which a triphosphate nucleoside transfers covalently its γ-phosphate to the enzyme, and the diphosphate nucleoside acceptor is subsequently bound to the same site to accept the phosphate.NDPK is known to phosphorylate a wide range of NAs bearing various structural modifications.Purified NDPK [27] is able to phosphorylate GDP (control) into GTP and also AT-8500 into AT-9010, the expected product (Fig 1B and Table 1).NDPK shows a high conversion rate and is only 4-fold slower with AT-8500 as substrate than GDP (1,300 ± 100 nmol/min/nmol protein versus 5,400 ± 800, respectively).We crystallized the NDPK with AT-9010 (Fig 3B).The crystal structure was solved at 1.9 Å resolution.A hexamer defines the asymmetric unit with 6 extra densities corresponding to 6 AT-8500 diphosphate.Interestingly, only the α and β-phosphates of the 5 0 -triphosphate moiety are visible.A new density is observed near the NE2 of His 118 of each monomer, which correspond to the histidine phosphorylated state His 118-P in Fig 3B (inset).The AT-9010 γ-phosphate is transferred to the latter due to the protein ability to catalyze the reaction in both directions via a ping-pong mechanism.Structural features of the nucleotide binding site are consistent with the required NDPK substrate promiscuity.Direct comparison with NDPK in complex with GDP (PDB: 1NUE) shows a conserved base-stacking between the guanine and Phe 60.The hydrogen bond between Lys 12 and the 3 0 -OH of the ribose is no longer present because of a displacement of the ribose, probably due to the 2 0 -C-methyl, 2 0 -F modification.This shift allows the 3 0 -OH to make direct contact with the phosphate group of His 118-P (3 Å).The β-phosphate is engaged in hydrogen bonds with Arg 88 and is facing the catalytic site of the enzyme.

Discussion
Here, we cast light on the activation pathway of bemnifosbuvir and its epimer AT-752.Our work quantitates the differential activation of these substrates relative to their natural counterpart GTP.We visualize and measure essential interactions these 2 compounds and their metabolites engage with specific enzymes involved in their sequential activation.
Turnover measurements (Table 1) show that the rate-limiting steps are generally occurring early in the pathway (CatA/CES1, HINT1, and ADALP1) rather than at the level of the kinases.On these slow, early steps, we note a de-esterification rate approximately 14-fold higher for AT-511 than AT-281 into the common intermediate AT-551 by CatA.Comparatively, CES1 reacts approximately 300-fold slower than CatA with AT-511 but shows an opposite stereopreference with AT-281 being preferred approximately 8-fold over AT-511.These results are aligned with those obtained with PSI-7851, a racemic mixture of PSI-7976 and PSI-7977 (SOF), which shows both the same CatA and stereoisomer preference [14].A number of NA prodrugs (e.g., SOF, TAF, RDV) need CatA/CES1 [10,11,37,38]-mediated activation.However, CatA is inhibited by the antiviral drug telaprevir (an HCV protease inhibitor) [10,14], an observation that casts light on drug-drug interaction at the molecular level.
Murakami and colleagues [14] have shown that HINT1 is the rate-limiting step along the SOF activation pathway.Its metabolic intermediate PSI-352707 (S1 Fig) binds poorly to HINT1, precluding saturating substrate conditions in infected cells and determination of kcat and Km in enzyme assays.In the case of AT-551, HINT1 is also the rate-limiting step (Table 1).Despite intense crystal soaking efforts using AT-551, we observed only weak positive partial density orthogonal to the expected, "natural substrate" position where phosphate of AT551 faces the catalytic His 112.Co-crystallization studies resulted in a clean HINT1:AT-8003 complex structure, which proposes a possible structural basis for the kinetic bottleneck along the SOF activation line [14].It may well be that the higher HINT1 conversion activity of AT-551 as compared to the corresponding intermediate metabolite of SOF, PSI-352707, which differ only by their nucleobase, results from the higher stacking power of purines versus uracil onto Ile 44.The HINT1 limiting step, however, does not seem to impact accumulation of the SOF 5 0 -triphosphate form (nor AT-9010) in cells.The compound pools are then readily available for selective use by the NS5B viral RdRp in HCV-infected liver cells.
These results indicate that cell types and their relative expression of these "early step" enzymes CatA/CES1 and HINT1 are expected to play a major role in drug activation: (i) the esterase-mediated deprotection step may well govern the overall rate of conversion of the drug into its active metabolite; (ii) the relative abundance of CatA-and CES1-like enzymes should influence these early steps; and (iii) further optimization in antiviral activity might be possible through optimization of the prodrug part in the 5 0 portion of the NA rather than focusing on optimizing the kinase steps.The choice of targeted infected cells may ultimately point to the most appropriate epimer and prodrug.In the case of ADALP1, unlike HINT1, the substrate and products are superimposable in the crystal structure, indicating that this complex should be reliable and useful for drug design.
Ideally, the NA activation line must be proven functional in cells, organs, or animal models supporting preclinical tests before commitment to full clinical trials in humans.AlphaFold2 [39]generated structural models of homologous animal enzymes did not point to obvious polymorphisms potentially deleterious for activation.The very high degree of structural conservation along species suggests a vital role in mammals for these activation enzymes.As such, monitoring differential expression of HINT1 should remain a sufficient and essential asset in predicting P-N hydrolysis in a given tissue or cell.The situation is somewhat similar for ADALP1, which shows significant polymorphisms across organisms, although mapping them invariably points to external, solvent exposed loops, with no direct impact on the NA binding site.
In the case of exposition of cells to bemnifosbuvir, its corresponding 5 0 -triphosphate AT-9010 has been demonstrated in cellular concentrations up to approximately 700 μM in human primary cells and cell types incubated in vitro with AT-511, including bronchial and nasal epithelial cells [18], hepatocytes, and Huh7 cells [24].Thus, none of the enzymes studied here limit the formation of the active compound AT-9010 in these cells.Bemnifosbuvir, though, shows reduced antiviral activity in some cell lines, e.g., in VeroE6, HeLa, or MRC-5 cells [18].We note that cynomolgus monkey hepatocytes have been reported to activate bemnifosbuvir to levels 50-fold lower than Huh7 cells [18].The superimposition of human and cynomolgus HINT1 and ADALP1 structures does not point to any polymorphism that could substantiate differences in substrate binding or catalysis.
Bemnifosbuvir shares 2 activation enzymes (HINT1 and NDPK) with SOF and RDV.Hence, the HINT1 structure should represent the most relevant enzyme to optimize the prodrug part of such antivirals.Two other enzymes ADALP1 and GUK1 are specific for diamino purine-and guanine-containing NAs, respectively.With the abacavir success in mind [22], the ADALP1 structure presented here should guide further exploration of chemical diversity at the N 6 -amino purine group.
A wealth of structural data already exists regarding NAs in complex with RNA virus replicases, e.g., HCV [15]; picornavirus [40]; Reovirus [41]; SARS-CoV-2 [13,28,42]; for a repository of stacked RdRp structures, see [43].Much like along the NA activation line, analysis of selectivity of viral polymerase in ternary complexes also identifies proximal amino acid side chains potentially able to accommodate chemical decorations without impacting the NA-5 0triphophate incorporation efficiency.
Finally, an NA drug potential is also conversely determined by its innocuity toward host enzymes that could misuse their 5 0 -triphosphate form and alter cell metabolism or host nucleic acid.Human mitochondrial RNA polymerase (POLRMT) inhibition accounts, in part, for cellular toxicity of NAs [44].Much like SOF-TP, AT-9010 is a poor substrate for POLRMT, which uses GTP >69,000-fold better [18].Comparative examination of ternary complexes of viral RdRps and POLRMT [45] with their RNA and NA-TP substrates may further document the strategy to rationally integrate NA drug design from chemistry to selective inhibition of viral growth [44].
Our work casts light on a whole set of previously ill-defined reactions involved in drug activation and contributes to the integration of NA drug-design from bioavailability to precise and selective mechanism of action as direct-acting antivirals.

Materials
Recombinant human CatA, cathepsin L (CatL), CES1, and trans-epoxy-succinyl-L-leucylamino(4-guanidino)butane (E-64) were purchased from R&D Systems.The respective genes coding for the human recombinant enzymes HINT1, ADALP1, GUK1, and NDPK were coded into a pNC-ET28 expression vector with an N-terminal hexahistidine tag (6xHis) and a TEV protease recognition site, purchased from Twist Bioscience.In humans, the NDPK gene is expressed as several isoforms, NDPK-A and NDPK-B being both in synthesis of nucleoside triphosphates other than ATP.The sequence of NDPK-B [26,27,46] was used throughout the study.Control compound AMP-NH 2 was from Biosynth, N 6 -Me-AMP was purchased from BioLog Life Science Institute, and AMP, IMP, GMP, GDP, and GTP were from Sigma Aldrich.Other compounds with a name starting with "AT" were synthesized by TopPharmAU : Pleasenotethata , USA.

Protein expression
The DNA sequence encoding for ADALP1, NDPK, GUK1, and HINT1 were synthesized and cloned into a pET-28a (+) vector for HINT1 and GUK1, into a pNC-ET28 vector for NDPK, and into a pET-28 for ADALP1 (Twist Biosciences).
All proteins were expressed in E. coli NEB C2566 cells (New England Biolabs).The cells were grown at 37˚C in TB medium for HINT1, GUK1, and NDPK and LB Broth for ADALP1 containing 50 μg/mL Kanamycin until the absorbance at 600 nm reached 0.6 to 0.8.At this stage, induction was started with 0.5 mM IPTG and the cells were grown overnight at 16˚C.Then, the cells were harvested and the pellets were suspended in lysis buffer (50 mM Tris (pH = 8.0), 300 mM NaCl, 10 mM imidazole, 1 mM PMSF, 0.25 mg/mL lysozyme, and 10 μg/ mL DNase) for HINT1, GUK1, and NDPK, before storing them at −80˚C.For ADALP1, the pellet was freeze-dried and stored at −80˚C.

HINT1, GUK1, and NDPK purification
The lysate was disrupted and cleared by centrifugation at 12,000g for 30 min at 10˚C.Then, the supernatant was loaded onto a Ni-NTA-beads (Thermo Fisher) batch and washed with 50 mM Tris (pH = 8.0), 300 mM NaCl, 20 mM imidazole buffer.The recombinant protein was then eluted with 50 mM Tris (pH = 8.0), 300 mM NaCl, 250 mM imidazole buffer.AboutAU : Pleaseno 10% glycerol was added for HINT1.TEV protease was used to remove the N-terminus His-tag during an overnight dialysis in a buffer (50 mM Tris (pH = 8.0), 150 mM NaCl, 1 mM DTT) at 4˚C for HINT1, at 37˚C for GUK1, and at room temperature for NDPK.
A second purification step on Ni-NTA beads was performed to remove any uncleaved protein before to achieve a SEC purification onto a Superdex 75 16/600 GE column (GE Healthcare) in a final buffer containing 10 mM Tris (pH = 8.0), 50 mM NaCl for HINT1; 20 mM Tris (pH = 8.0), 100 mM NaCl, 1 mM DTT for GUK1; and 50 mM Tris (pH = 8.0), 150 mM NaCl, 1 mM DTT for NDPK.The fractions were analyzed on SDS-PAGE and the ones containing the pure target protein were pooled.The purified protein was then concentrated until 10 mg/ mL, aliquoted, and stored at −80˚C.Same protocol was used to express and purify uncleaved NDPK for activity assays.
TEV protease was used to remove the N-terminus His-tag during an overnight dialysis in a buffer 50 mM HEPES (pH = 7.5), 300 mM NaCl, 10% glycerol, 1 mM DTT.The untagged ADALP1 was further purified by a Talon Superflow beads batch (Merck).
Finally, ADALP1 was purified by size exclusion chromatography using a S75 16/60 GE Column (GE Healthcare) equilibrated in buffer containing 10 mM Tris (pH = 8.0), 40 mM NaCl, 1 mM TCEP.The single ADALP1 peak corresponded to the monomeric form.It was concentrated to 10 mg/mL, flash-frozen in liquid nitrogen, and stored at −80˚C.

CatA enzyme assay
Human recombinant CatA was first activated by following the manufacturer's instructions (R&D Systems), incubating at 37˚C 10 μg/mL CatA with 1 μg/mL CatL in an activation buffer (25 mM MES (pH = 6.0), 5 mM DTT) during 30 min.AboutAU : PleasenotethatasperPLOSstyle; numer 10 μM E-64 (CatL inhibitor) was then added before aliquoting and storage of activated CatA at −80˚C.AT-511 and AT-281 hydrolysis by activated CatA was measured by incubating 100 μM compound in a reaction buffer containing 25 mM MES (pH 6.5), 100 mM NaCl, 1 mM DTT, 0.1% NP-40, and 20 nM enzyme at 37˚C for 45 min.The reaction was started by adding the enzyme.At various time points, 10 μL aliquots were collected from the reaction mixtures mixed with EDTA 10 mM (final concentration) and stopped by heating the samples at 95˚C for 5 min.The sample were filtered on centrifugal filters Microcon-10 (Sigma Aldrich).The filtrates were mixed with triethylammonium bicarbonate (TEAB) 1M (1:1) and analyzed onto a C18 reverse phase column (2.5 μm, 4.6 by 100 mm, Xbridge C18 BEH, Waters) equipped with a guard column; substrates and products were detected at 260 nm and eluted at a flow rate of 1 ml/min at 35˚C with a gradient of 0% to 4% of MeCN in TEAB 50 mM (pH = 7) in 10 min.

CES1 enzyme assay
AT-511 and AT-281 hydrolysis by human recombinant CES1 (100 nM) was assayed by incubating at 37˚C the enzyme with 100 μM compound in 50 mM Tris buffer (pH = 7.5), 0.1% NP-40 and 1 mM DTT for 2 h.The reaction was started by adding the substrate.At various time points, 20 μL aliquots were collected from the reaction mixtures mixed with EDTA 10 mM (final concentration) and stopped by heating the samples at 95˚C for 5 min.The sample were filtered on centrifugal filters Microcon-10 (Sigma-Aldrich).The filtrates were mixed with TEAB 1M (1:1) and analyzed onto a C18 reverse phase column (2.5 μm, 4.6 by 100 mm, Xbridge C18 BEH, Waters) equipped with a guard column; substrates and products were detected at 260 nm and eluted at a flow rate of 1 ml/min at 35˚C with a gradient of 0% to 4% of MeCN in TEAB 50 mM (pH = 7) in 10 min.

HINT1 enzyme assay
The hydrolytic reactions with HINT1 (100 nM) were performed on AMP-NH 2 and AT-551 (200 μM) in 20 mM HEPES buffer (pH 7.2), 20 mM KCl, 1 mM MgCl 2 , and 1 mM DTT, at 37˚C for 2 h.The reaction was started by adding the enzyme.At various time points, 10 μL aliquots were collected from the reaction mixtures mixed with EDTA 10 mM (final concentration) and stopped by heating the samples at 95˚C for 5 min.The sample are filtered on AcroPrep Advance 96-Well Filter Plates with 3K Omega membrane (Pall).The filtrates were mixed with TEAB 1M and analyzed onto a C18 reverse phase column (3 μm, 3 by 150 mm, Acclaim Polar Advantage II, Thermo Fisher) equipped with a guard column; substrates and products were detected at 260 nm and eluted at a flow rate of 0.5 ml/min at 30˚C with a gradient of 0% to 15% of MeCN in TEAB 50 mM (pH = 7) in 17 min.

ADALP1 enzyme assay
ADALP1 (100 nM) activity was measured on several substrates (N 6 -Me-AMP, AT-8003, AT-8002, AT-8004, AT-8010, AT-551, AT-259, AT-229) at 200 μM concentration in a reaction buffer containing BTP (pH = 6.8), 100 mM NaCl, and 1 mM DTT.Reactions were incubated at 37˚C for 2 h.AboutAU : PleasenotethatasperPLOSstyle; numeralsarenotallowedatthebeginningofasen 1 μM ADALP1 was also tested with AT-8002, AT-551, AT-229, and AT-259.The reaction was started by adding the substrate.At various time points, 10 μL aliquots were collected from the reaction mixtures diluted with H 2 O and stopped by heating the samples at 95˚C for 5 min.The sample were filtered on AcroPrep Advance 96-Well Filter Plates with 3K Omega membrane (Pall).The filtrates were mixed with TEAB 1M and analyzed onto a C18 reverse phase column (2.5 μm, 4.6 by 100 mm, Xbridge C18 BEH Premier, Waters) equipped with a guard column; N 6 -Me-AMP, AT-8002, and products were detected at 260 nm and eluted at a flow rate of 1 ml/min at 35˚C with a gradient of 0% to 10% of MeCN in TEAB 50 mM (pH = 7) in 10 min.AT-8003, AT-551, and products were eluted at a flow rate of 1 ml/ min at 35˚C with a gradient of 0% to 15% of MeCN in TEAB 50 mM (pH = 7) in 10 min.AT-8004, AT-8010, AT-229, AT-259, and products were eluted at a flow rate of 1 ml/min at 35˚C with a gradient of 0% to 25% of MeCN in TEAB 50 mM (pH = 7) in 10 min.

GUK1 enzyme assay
GUK1 (40 nM) was assayed on GMP, AT-8003, and AT-8001 (200 μM) in a reaction buffer containing 50 mM Tris buffer (pH = 8.0), 50 mM KCl, 5 mM MgCl 2 , 1 mM ATP, and 1 mM DTT.Reactions were incubated at 37˚C for 1 h.The reaction was started by adding the enzyme.At various time points, 10 μL aliquots were collected from the reaction mixtures, mixed with EDTA 10 mM (final concentration), and stopped by heating the samples at 95˚C for 5 min.The sample were filtered on AcroPrep Advance 96-Well Filter Plates with 3K Omega membrane (Pall).The filtrates were mixed with TEAB 1M and analyzed onto a C18 reverse phase column (3 μm, 3 by 150 mm, Acclaim Polar Advantage II, Thermo Fisher) equipped with a guard column and equilibrated with TEAB 50 mM (pH = 7) (buffer A).Substrates and products were detected at 260 nm and eluted at a flow rate of 0.5 ml/min at 30˚C with a gradient of 0% to 10% of MeCN in TEAB 50 mM (pH = 7) in 12 min.

NDPK enzyme assay
NDPK (20 nM) was assayed on GDP and AT-8500 (200 μM) in a reaction buffer containing 50 mM Tris buffer (pH = 8.0), 50 mM KCl, 5 mM MgCl 2 , 1 mM ATP, and 1 mM DTT.Reactions were incubated at 37˚C for 1 h.The reaction was started by adding the enzyme.At various time points, 10 μL aliquots were collected from the reaction mixtures mixed with EDTA 10 mM and stopped by heating the samples at 95˚C for 5 min.The samples were filtered on Acro-Prep Advance 96-Well Filter Plates with 3K Omega membrane (Pall).The filtrates were mixed with TEAB 1M and analyzed onto a C18 reverse phase column (3 μm, 3 by 150 mm, Acclaim Polar Advantage II, Thermo Fisher) equipped with a guard column; substrates and products were detected at 260 nm and eluted at a flow rate of 0.5 ml/min at 30˚C with a gradient of 0% to 10% of MeCN in TEAB 50 mM (pH = 7) in 12 min.

Crystallization of HINT1
Crystallization conditions were adapted from [21].Briefly, co-crystals with AT-8003 were grown at 293.15 K, using a 1:1 ratio of HINT1 at 10 mg/mL with AT-8003 (final concentration 25 mM) to precipitant solution 0.1 M MES (pH 6.1 to 6.5), 27% to 30% PEG 8000.Crystals grew in a few days and were cryo-protected with reservoir solution supplemented with 20% PEG 400, and flash-frozen in liquid nitrogen at 100 K.

Crystallization of ADALP1
Crystallization assay were set up using the sitting drop vapor diffusion method, using a 2:1 ratio of protein solution to the reservoir solution at 293.15 K.The co-crystals with AT-8001 were obtained by mixing ADALP1 at 17 mg/mL with AT-8003 (final concentration 10 mM) and 1:10,000 w/w dilution a-Chymotrypsin prior to crystallization and were grown from 1.1 to 2.1 M ammonium sulfate, 0.1 M sodium cacodylate/HCl (pH 5.3 to 6.3), and 0.2 M sodium chloride for 2 moAU : Pleasenotethatallunitsoftimeði:e:; month; day; etc:Þhavebeenabbreviatedtoenforce .Crystals were cryo-protected with reservoir solution supplemented with 20% glycerol and flash-frozen in liquid nitrogen at 100 K.

Crystallization of GUK1
Large co-crystals of GUK1 (17 mg/mL) with AT-8001 (2 mM) supplemented with 5 mM MgCl2 were grown using the sitting drop method with a 3:1 ratio of protein solution to precipitant solution consisting of 66 mM sodium cacodylate (pH 6.5 to 7.5), 12.2% to 22.2% PEG 3350, 8.3% PEG 4000, 33 mM MES (pH 6.5), and 66 mM magnesium chloride over a period of 10 d at 293.15 K. Crystals were cryo-protected with reservoir solution supplemented with 20% glycerol and flash-frozen in liquid nitrogen at 100 K.

Crystallization of NDPK
Crystallization conditions were adapted from [46] using the sitting-drop vapor diffusion method using crystallization buffer.All crystals were grown at 277.15 K, using a 1:1 ratio of protein (5 mg/mL) mixed with AT-9010 (final concentration 1 mM) to precipitant solution (12% PEG 4000, 50 mM Tris (pH = 8.4), 16% glycerol, 1 mM DTT).Crystals grew in 3 d and were flash-frozen directly in liquid nitrogen at 100 K.

Structure determination of the human HINT1:AT-8003 complex
The dataset of HINT1 in complex with AT-8003 was collected on the Proxima-1 beamline at the Synchrotron SOLEIL.Dataset was processed using AUTOPROC [47].The phase was obtained using Molecular Replacement using MOLREP [48] with the PDB entry 6N3V as a model.The ligand AT-8003 and geometry description files were generated using the GRADE2 server [49].Structure handling and refinement were done using COOT [50] and BUSTER [51].

Structure determination of the human ADALP 1:AT-8001 complex
The dataset of ADALP 1 in complex with AT-8001 was collected on the Proxima-2 beamline at the Synchrotron SOLEIL.Dataset was processed using AUTOPROC [47], and processed data were sent to the CCP4 cloud suite [52].The pipeline is fully automatic.The phase was obtained using Molecular Replacement-MORDA [52] using PDB entry 6IV5 as a model.The ligand AT-8001 and geometry description files were generated using the GRADE2 server [49].Minor correction on the structure and further refinement were done using COOT [50] and REFMAC5 [53], respectively.

Structure determination of the human GUK 1:AT-8001 complex
The dataset of GUK1 in complex with AT-8001 was collected on the ESRF ID23-1 beamline.Dataset was processed using AUTOPROC [47].The phase was obtained using Molecular Replacement-MORDA [52] using PDB entry 4F4J as a search model and optimized with ARP/WARP [54].Structure refinement was done using COOT [50] and BUSTER [51], respectively.

Structure determination of the human NDPK:AT-8500 complex
The dataset of NDPK in complex with AT-8500 was collected on the Proxima-2 beamline at the Synchrotron SOLEIL.Dataset was handled using the CCP4 suite [55].Images were processed by XDS [56] and AIMLESS [57].The phase was obtained using Molecular Replacement-PHASER [58] with the PDB entry 1NUE as a model.The ligand AT-8500 corresponding to the diphosphate form of AT-9010 was generated using the AceDRG program [59].Structure handling and refinement were done using COOT [50] and REFMAC5 [53] software, respectively.
All electron-density maps were inspected using COOT [50].Extra density accounting for ions, and/or compounds were observed for all complexed structures.The structures were evaluated using MOLPROBITY [60] and PROCHECK [61].Structural analysis and high-resolution figures were done with UCSF ChimeraX [62].Facilities, statistics of data collections, refinements, and PDB deposition code are given in Table 2.

Modeling of GUK 1 closed conformation
Sequence of human GUK 1 was submitted to an homology modeling process using modeler V10 and the structure of mouse GUK complexed with GMP and ADP (PDB 1LVG) as template.The 2 sequences have 88% identity, which makes us confident that the resulting closed conformation model is reasonable.

Docking of AT-8003 onto HINT1
The 3D model of HINT1 was energy minimized by the steepest gradient method of energy minimization followed by conjugate gradient minimization, using the MMTK and AMBERpackages [63][64][65].Mol2 and PDB files format of the ligands and receptor were converted to PDBQT format using CHIMERA prior to docking.All the water and solvent atoms of the protein were removed, and the polar hydrogen and polar charge were added onto the ions and ligand prior to docking.The protein was kept rigid, while the ligand was allowed to rotate and explore more flexible binding pockets.Docking of the respective ligands into the cavity were performed iteratively using AUTODOCK VINA-version 1.1.2[66].The best poses from the first round of docking were used as seed for the second round.The resulting first round of docking were carefully analyzed to retain the best poses.The grid box size dimensions was designed to include the long binding cleft and the catalytic site; its dimensions were 37.87 × 16.12 × 20.89.The default scoring function was used for docking.
Ten binding modes of the docked complexes were obtained and sorted based on their binding energy, and amino acid residues present at a distance less of 3 Å were considered as the binding partners of the ligands.Six binding modes with the phosphate oriented toward the catalytic site were kept and compared to the experimental structure.One binding mode is overlapped with the structure.A control experiment with AMP was also performed.Figures representing the docked complexes have been generated using CHIMERA [62].

Docking of GMP and AT-8001 in GUK1 closed conformation
The 3D model of GUK 1 was energy minimized by the steepest gradient method of energy minimization followed by conjugate gradient minimization, using the MMTK and AMBER packages [63][64][65].Mol2 and PDB files format of the ligands and receptor were converted to PDBQT format using chimera prior to docking.All the water and solvent atoms of the protein were removed, and the polar hydrogen and polar charge were added onto the ions and ligand prior to docking.The protein was kept rigid, while the ligand was allowed to rotate and explore more flexible binding pockets.Docking of the respective ligands into the cavity were performed iteratively using AUTODOCK VINA [66].The best poses from the first round of docking were used as seed for the second round.The resulting first round of docking were carefully analyzed to retain the best poses.The grid box size dimensions were first 40 × 40 × 40, to verify that our ligands will preferentially bind in the catalytic site.The grid box size was further optimized to 23.2 × 25.6 × 21.2, thus covering the binding pockets; the default scoring function was used for docking.As control for the procedure, GMP was docked following the same protocol, and the final pose is virtually identical to the one measure in the experimental structure (PDB 1LVG).
Binding modes of the docked complexes were obtained and sorted based on their binding energy; ions and amino acid residues present at a distance less of 3 Å were considered as the binding partners of the ligands.Binding modes were compared to theses of the native structure.The interaction figures representing the docked complexes have been generated using CHIMERAX [67].

Fig 2 .
Fig 2. (A) Structure of human HINT1 in complex with AT-8003.Homodimer of human HINT1 structure complexed with 1 AT-8003 compound is represented in ribbon and transparent surface, while compound is represented in sticks with heteroatom colors.On the right side, an enlargement of the molecule in the active site is presented at 90˚.Side chains of interacting residues with the molecule are shown in sticks.The F o -F c omit map of the molecule is presented at 1 σ.(B) Binding mode involving direct coordination of the 2 0 and 3 0 O ribose of AMP (light green-PDB: 5ED6) by Asp 43. (C) Comparison of the relative

Fig 3 .
Fig 3. (A) Structure of human GUK1 in complex with AT-8001.Human GUK1 structure is represented in ribbon and transparent surface, while compound is represented in sticks with heteroatom colors.The GUK1 structure corresponds to the open conformation, and molecule is binding to the loading site.On the right side, a zoomed part of the molecule is presented in the loading site.Water and side chains of interacting residues with the molecule are shown in sticks.The F o -F c omit map of the molecule is presented at 1 σ.(B) Structure of human NDPK phosphorylated on H118 and in complex with AT-8500.Human NDPK-B structure (see Materials and methods) is represented in ribbon and transparent surface, while compound is represented in sticks with heteroatom colors.On the right side, a zoomed part of the molecule is presented in the active site in 2 orientations differing by 90˚.Side chains of interacting residues with the molecule are shown in sticks.The F o -F c omit map of the molecule as well as the transferred phosphate onto the catalytic histidine is presented at 1 σ.The data underlying this figure can be found at https://doi.org/10.2210/pdb8pts/pdband https://doi.org/10.2210/pdb8pie/pdb.https://doi.org/10.1371/journal.pbio.3002743.g003

considering 1
site model.The data underlying this supporting figure can be found at https:// zenodo.org/records/12606239.(TIF) S9 Fig. Docking studies of ligands onto HINT1.(A, B) Two views (side and tilted) visualizing the criteria for docking experiment: HINT1 structure is presented as the accessible surface (colored wheat).The green box represents the docking grid encompassing the binding and catalytic site; AT-8003 in predocking position (colored in blue) is outside the computational grid and far from the protein target.(C) Docking result overlapping all 10 poses.Docked molecules are colored in light green.None of the molecules are binding following the binding cleft.(D) Similar to C except only remains the 6 poses for which the phosphate is oriented toward the catalytic site.(E) The one pose that is similar to the experimental data.(F, G, H) Superimposition of the docked AT-8003 with the AT-8003 (colored in light pink) measured by crystallography.Phosphate, ribose, and part of the base are fully superimposed, and only the extremity of the base is shifted.(I) Control of docked AMP (colored in purple) showing that the molecule is fully positioned within the binding site.(TIF) S1 Raw Images.Raw blots for S2 Fig. (PDF)

Table 2 . Data-collection and refinement statistics for native human enzymes of the bemnifosbuvir activation pathway complexed with substrate or product compounds.
*Outliers are in highly flexible regions of the protein.https://doi.org/10.1371/journal.pbio.3002743.t002